Analytical calculation of the drag force near drag crisis of a falling sphere 



Armando V.D.B. Assis, M.H.R. Tragtenberg, and N.S. Branco 
Departamento de Fisica, 
Universidade Federal de Santa Catarina, 
88040-900, Florianopolis, SC, Brazil 

(Dated: September 9, 2010) 

We obtain analitically the v 2 dependence of the drag force on a falling sphere close to the drag 
crisis, as well as the drag coefficient at the drag crisis, with excellent agreement with experiment. 
We take into account the effects of viscosity in creating a turbulent boundary layer and perform the 
calculations using the Navier-Stokes equation. 
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We report in this paper the first derivation of both 
the v 2 dependence of the drag force and the drag coeffi- 
cient at the drag crisis (abrupt decrease of the drag co- 
efficient) for a falling sphere, which is in excellent agree- 
ment with experimental results and is derived from the 
Navier-Stokes equation. Previous works are far from the 
well-known experimental results [5] . We remark that 
this is, to our knowledge, the first analytical derivation 
taking into account the Navier-Stokes equation and the 
character of the boundary layer, and which is in excellent 
agreement with experimental data. 

There is no theoretical determination of the velocity 
dependence of the drag force based on the Navier-Stokes 
equation that take into account the physics of the bound- 
ary layer and accomplish with the dependence of the drag 
force on the sphere velocity close to the drag crisis. We 
understand that these achievements are important steps 
towards the full understanding of the underlying aspects 
of the physics of the drag crisis. The works in the liter- 
ature that support our results are [Tl-|13|. in agreement 
with the result for the drag coefficient reported in this 
paper. 

We start considering the case of the sphere falling in 
the stationary state, with constant velocity in relation 
to the ground. We consider the reference frame fixed 
on the ground, where the Navier-Stokes equation for an 
isotropic fluid reads 
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where v is the fluid velocity field, p is the fluid density, g 
is the external gravitational field, m is the sphere mass, 
p is the pressure field, rj and £ are the first and the sec- 
ond viscosity coefficients taken as constants and the dot 
denotes the total time derivative. The Einstein summa- 
tion convention is being used on repeated indices. For 
divergence free velocity field fluids (constant density is a 
sufficient condition), the Navier-Stokes equations turns 



pv — pg + Vp — ri\7 v = 0. 
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When we go to the reference frame fixed on the sphere, 
the Navier-Stokes has the same form, but now g — > 
—F/m, where F is the fluid force on the sphere at rest. 
We assume null fluid velocity at the sphere surface, the 
nonslip boundary condition on the surface of the sphere, 
and steady state fluid velocity —h(t)e z far away from the 
sphere, where h(t)e z is the the fall velocity of the cen- 
ter of the sphere in relation to the ground and e 2 is the 
vertical upwards versor. 

We want to determine the force F the fluid exerts on 
the sphere. Back to the general case, from the continuity 
and Navier-Stokes equations we arrive to the equation 
below for each component: 
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Integrating Eq.4 in a fixed, arbitrary and nondeformable 
control volume, we have: 
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where n is the tensor: 
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1 is the identity tensor and v(£>vis the dyadic product. In 
general case, —F/m is the non-inertial acceleration field, 
the same at any point of the non-inertial frame fixed on 
the sphere at a given instant. Hence, the fluid force on 
the sphere is given by 



m Jcv P dv V/cs at J cv 
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where CV and CS are the control volume and the control 
surface, respectively. 

The fluid velocity field may be divided in three main re- 
gions: the mainstream, the boundary layer and the wake 
in Fig. 1. 
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Figure 1: Schematic draw of the main regions. 

The stationary velocity field in the mainstream Voo(f) is 
taken as irrotational in virtue of the rapid decrease of 
deformation outside the boundary layer compared to the 
internal boundary layer region, at high Reynolds num- 
bers. Hence, the velocity field is taken as potential flow 
Vooi'r) = Vif(r) outside the boundary layer. For an in- 
compressible fluid, K{r) must obey the Laplace equation. 
The boundary conditions for the stationary velocity field 
are: at infinity it is lim Voo(f) — —h°°(t)e z , whereas 

|r|— >oo 

it is tangent to the outer surface of the boundary layer. 
The solution for the stationary fluid velocity field at the 
mainstrean is 



«oo(r) = -h°°(t) 
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where R is the sphere radius and 6 is the thickness of 
the boundary layer. The mainstream stationary velocity 
field obeys the Navier-Stokes equation in the form of Eq. 
(3) (with g — > —F/m). Since V x F = 0, we define a 
potential <p°°. Thus, we find a Bernoulli field, such that: 



(9) 



where is the pressure scalar field in the mainstream 
region, pjJo a constant pressure field. In order to evaluate 
the force F°° the fluid exerts on the sphere we recall the 
Eqs. (6), (7), in the stationary state 
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We need to define an appropriate control volume and its 
associated control surface. We choose the control surface 
as the AFGBA surface, where BA is in contact with the 
wake, as can be seen in Fig. 2. 



e-e d 




(10) 



Figure 2: Schematic draw of the control region. 

The time dependent term in the rigth hand side of Eq. 
(10) can be determined under the assumption that we 
are interested in the brink of the drag crisis, where the 
boundary layer is in a turbulent transition regime. In 
that regime, the turbulent boundary layer does not have 
a stationary state. But we are interested in the average 
effect, which can be taken into account by considering 
the ensemble average, which is time independent. Then, 
we can neglect this term. 

We can also neglect the term of Eq. (11), con- 
sidering that the turbulent boundary layer has a bigger 
velocity gradient in the surface AB (AB in contact with 
the boundary layer) in comparison with the velocity gra- 
dient in FG (the onset profile of the velocity internal to 
the boundary layer will be, in few lines below, modeled 
by a step function). This is because we are close to the 
drag crisis, and the nonslip boundary condition imposes 
big changes in velocity from the surface to regions close to 
the sphere, within the boundary layer, and we are able to 
choose the surface FG as far as we want from the sphere 
surface, in order to have a smaller velocity gradient, even 
though this gradient is within a turbulent region. 

The integral in the denominator of Eq. (10) is m + 
rriBLi t ne sum of the mass of the sphere with the mass 
of the boundary layer. The tensor Vaa ® Voo is null at 
the surface AB in contact with the wake, by the non- 
slip boundary condition, as well as at the surface close to 
the separation points (AF and GB in Fig. 2), since the 
velocity is close to zero in this region. We end up with 
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In order to calculate we have to determine the 

steady state pressure field and the tensor (Uoo ® «oo)- 
Concerning the pressure field for a laminar boundary 
layer, we have no change in this field for different dis- 
tances from the surface of the sphere 
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PP. However, for a turbulent boundary layer, close to the 
drag crisis, Eq. (13) is not true for the pressure field p, 
but it is a good approximation for the pressure field p — 
p+pipp/m. This pressure field incorporates the potential 
related to the non-inertial irrotational force field. We 
can assume this in the region between the surfaces FG 
and CD, since this region is sufficiently thin and laminar. 
Then, 
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The pressure field in the GBAF surface (BA in con- 
tact with the wake) will be constant [31 [5] . Then, if we 
are able to calculate the pressure field at the separation 
point S (= B), we can calculate it in that surface. We 
invoke that we can draw an SG (= BG) current line, at 
the drag crisis, close to the B and G points where we 
have low shear, since it is the separation point neighbor- 
hood. Throughout this line, Euler equation is valid and 
Bernoulli law too. Then we are able to support that 
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since 8 << R, where S is the radial distance from the 
sphere surface until the FG surface. Then the pressure 
field close to wake and the separation point, is given by 
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We will model the square velocity profile inside the 
boundary layer as a step function, since it is a high gra- 
dient velocity field close to the drag crisis, and we are able 
to expand the time average of this field into a Fourier se- 
ries. The value of the square velocity at the distance 5' 
is then given by 



where a is the angle between the velocity of the fluid 
and the normal vector of the FG spherical surface 
(see Fig. 2). Then, the tensor which gives the lin- 
ear momentum flux through FG is (v (§5 v) FG ■ h — 
Vfg (R + S' ,9,t)[vpG (R + 5' ,6,t) ■ n], and has an en- 
semble average given by 
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Then, the pressure field close to the separation point will 
be 
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where GBAF is the region composed by the wake in con- 
tact with the sphere plus the neighboring region to the 
separation point. Since the region FG is turbulent, the 
velocity field is not tangential there. We have to perform 
a time average there. The turbulent velocity field, there, 
is 

vfg — vfg (R + S', 9, t) [cosa(t) e r + sina(£) eg] , (19) 



(20) 

Assuming randomness in the angles a between the fluid 
elements and the normal of the sphere surface, we find 
that 



(ilao ® Woo) ■ n = - (v 2 FG (R + S', 9, t)) t e r 
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From the Eqs. (14), (18), (21) and (12) we are able to 
calculate the force the fluid exerts on the sphere 
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The drag force is eventually given by 
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where is the viscous drag force, since F°° /m = —g in 
the stationary state. The drag coefficient Cd at the drag 
crisis is 



Cn = 



2F a 



= 0.44, 
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since 9 s = 70,4° is the separation angle [3]. This coeffi- 
cient is in excellent agreement with experimental results 
IB- 



ACKNOWLEDGMENTS 

The authors would like to thank CNPq and CAPES 
for partial financial support. 



2nd ed., english Edition. 536 pp. 

[2] J. E. Goff, Physics Education 39, 496 (2004). 

[3] H. Schlichting, Boundary-Layer Theory (McGraw-Hill 
Books, 1979), 2nd ed., english Edition. 

[4] O. Ladyzhenskaya, The mathematical theory of viscous 
imcompressible flow. Mathematics and its Appl, second 
volume (Gordon and Breach, New York-London-Paris, 
1969), 2nd ed., english Edition. 

[5] M. J. Chorin A., A Mathematical Introduction to Fluid 
Mechanics (Springer, 1993), 1st ed. 

[6] K. B. Arnold V.I., Topological methods in hydrodynam- 
ics, Applied Mathematical Sciences (Springer- Verlag New 
York, 1998). 

[7] S. W. V. S. Michael R Smith, David K. Hilton, Physics 

of Fluids Letters 11, 751 (1999). 
[8] S. W. V. S. Yeon Suk Choi, Michael R. Smith, National 

High Magnetic Field Laboratory 571, 66 (2001). 
[9] K. S. George Constantinescu, Physics of Fluids 16, 1449 

(2004). 

[10] W. Yeung, Journal of Wind Engineering and Industrial 

Aerodynamics 95, 253 (2007). 
[11] H. C. Jin Choi, Woo-Pyung Jeon, Physics of Fluids 18, 

041702/01 (2006). 
[12] R.-H. Y. M.-C. W. A.-B. W. Ming-Hsun Wu, Chih- 

Yung Wen, Journal of Fluid Mechanics 515, 233 (2004). 
[13] E. Achenbach, Journal of Fluid Mechanics 54, 465 

(1972). 



[1] L. Landau and E. Lifshitz, Course of Theoretical Physics, 
Vol.6. (Pergamon Press, London-Paris- Frankfurt, 1959), 



